A self-consistent perturbative evaluation of ground state energies: application to 

cohesive energies of spin lattices 
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The work presents a simple formalism which proposes an estimate of the ground state energy from 
a single reference function. It is based on a perturbative expansion but leads to non linear coupled 
equations. It can be viewed as well as a modified coupled cluster formulation. Applied to a series 
of spin lattices governed by model Hamiltonians the method leads to simple analytic solutions. The 
so-calculated cohesive energies are surprisingly accurate. Two examples illustrate its applicability 
to locate phase transition. 



I. INTRODUCTION 

The quantum many-body problem has been perfectly 
clarified a long time ago for the single-reference ap- 
proaches. Its logics is given by the linked cluster 
theorem ; 1 ! 2 ^ which relies on a perturbative expansion 
and its diagrammatic representation. A panoply of com- 
putational methods satisfying the basic size-consistency 
requirements are available. One of them is the Moller- 
Plesset order-by-order perturbative expansion^ which 
suffers from its slow convergence (and in some cases high- 
order divergences^) . Other methods lead to non-linear 
equations coupling the coefficients of the vectors which 
are singly or doubly excited with respect to the single 
reference <f>o. Among them one may quote the Cou- 
pled Electron Pair Approximation (CEPA) formalism§iL& 
which eliminates unlinked diagrams from a Singles and 
Doubles Configuration Interaction (SDCI) and intro- 
duces part of (in CEPA-3^) or all (in the SC 2 SDCI 
version^) of the Exclusion Principle Violating diagrams. 
The Coupled Cluster Method (CC M'li 11 . 12 ) has a better 
standard since it is an elegant and generalizable formal- 
ism. This method gives an exponential structure to the 
wave operator fl which transforms a reference function 
<&o in the exact solution $fo, = ^$0, ^ = exp(S). In 
practice S in limited to a certain set of excitations. The 
simplest version, called CCSD, introduces single and dou- 
ble excitations on the top of <&o and takes into account 
the connected effect of the quadruples. It does notes treat 
correctly the effect of the Triples, which has to be incor- 
porated perturbatively in the CCSD(T) formalism piii4 
The present work returns to a low-order semi- 
perturbative development from <j> - It leads to a system 
of equations fixing the coefficients of the singly and dou- 
bly excited vectors, which 

- does not require to assume an exponential structure 
of the wave-operator, 

- introduces the effect of deviations from additivity 
of excitation energies (deviations which is implicitly 
supposed to be negligible in CCSD), 

- and introduces some high-order effects through 
EPV corrections to the excitation energies. 



The method is applied to a few infinite periodic spin lat- 
tices governed by Heisenberg Hamiltonians which only 
consider interactions between nearest neighbors. In these 
problems if all bonds are identical, the solutions of our 
self-consistent perturbative (SCP) method lead to sim- 
ple polynomial equations. The results, compared to the 
exact solutions, are of surprising quality. The method is 
then employed to identify phase transitions in two differ- 
ent 2-D problems. 



II. FORMALISM 

A. Generalities 

Let us call <f>o a suitable single reference supposed to 
be a relevant approximation to the ground state of the 
system governed by an Hamiltonian H and $ m vectors 
orthogonal to <f>o • In the intermediate normalization 



#o = $o + Yl Cm ®r. 



(1) 



the exact energy is given by the eigenequation relative to 
$ , ($q\H - E\* ) =0 



E = Hi 



oo 



(2) 



where ffy = (pi\H\$j). The Hamiltonians is at most bi- 
electronic. One may call hi and hi the mono-electronic 
and bi-electronic parts of the Hamiltonian H = hi + hi. 
One may label m,n... (holes) the mono-electronic func- 
tion occupied in <j> and r,s... (particles) those which are 
not occupied in $o- 

In the following we shall refer to an Epstein-Nesbefei^iSiii 
zero-order Hamiltonian which is the diagonal part of 
the Hamiltonian in the basis of the N-electronic vectors 
{$ ,--$m--} 



The perturbation operator V 

V = H -H a 



(3) 



(4) 
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is zero-diagonal in this basis 

y = EE'i^>^i^>^i- ( 5 ) 

i 3 

The Hamiltonian being at most bi-electronic, it appears 
from Eq. 2 that the knowledge of the coefficient of the 
singly and doubly excited vectors which interact with <£> 
is sufficient to give the ground state energy. We shall call 
first generation S\ the set of vectors $i interacting with 

Si = {**}, (* < |ff|$ >? fe Q 

and T+ the operators creating the vectors <&j from $ , 

= T+§ . 

These operators may be written as 

I^ = o+a (l-*((rHa») (6) 

or 

T+ = a+a+a b a a (l - 6((rs\h 2 \ab))) (7) 

where S is the Kroneker symbol. The coefficients of 
the first generation vectors can be estimated from the 
eigenequations relative to them. For $j the equation 
($i\H — E\ty ) = can be written using Eq. 2 as 

+ J2 + E H *» c " = °- ( 8 ) 

j&Si a^Si 

To determine the coefficients C; of the first generation 
vectors, it is sufficient to have an estimate of the coeffi- 
cient of vectors $ Q of the second generation, i.e., those 
which interact with the vectors of Si. Among the vectors 
$ a belonging to the second generation one may distin- 
guish those which are obtained from $j by the actions 
of operators possible on <f>o (type 1), and the others 
which are obtained from $j by operators different 
from the T'+ (type 2). 

If one calls T the sum of the operators and of their 
adjoints 

r=E( T fc ++T * +x ) ( 9 ) 

k 

the operators R 

R=i-j2mm-T, (io) 
i 

is the sum of all the operators changing one or two mono- 
electronic functions which are different from T£ opera- 
tors. Notice that (R+$ \H\$ ) = 0. With the Hamilto- 
nians considered hereafter the second generation vectors 



$ a are either of type 1 or of type 2, i.e., are generated 
either as 

$ Q = T+T+* (type 1) 

or as 

$ Q = R+T+$ (type 2) 
as pictured in Fig. 1. 



(type 1) (type 2) 



Second generation 
First generation 
Reference 




FIG. 1: Illustration of the genealogic generation of from 
the reference function <E>o 



B. Evaluation of the coefficients of second 
generation vectors 

Regarding the vectors of type 1, the operators TJ7" are 
in general possible on but some of them are impos- 
sible (T+-5; = 0). Let call EPV(i) (Exclusion Principle 
Violating correction) the quantity^ 

EPV(i)= J2 H okCk- (11) 

k 

Noting that if <£> Q = Tj^^%, and if Tj" is a double excita- 
tion Hi a = -Hofc- Introducing Eq. 11 the Eq. 8 can be 
written 

{Ha - H m - EPV{i)) Ci +H M +J2 
+^ H k ^C T + $ . — CiCkj 

k 

+ y, @mR+^)c Rt ^ =o. (12) 

rn 

A first order evaluation of d is of course the Epstein- 
Nesbet 15 i 16 i 17 one 

Ci = Hio/(Hoo — Ha). (13) 

Introducing EPV corrections in the energy denominator 

d = H i0 /(H 00 - H u + EPV(i)), (14) 
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incorporates an infinite summation of diagrams One 
may use a second order perturbation theory to evaluate 



the coefficients C, 



and C 



The vector of type 1 T+$; = T+$ fe = T+T+<5> can also 
be reached from <f>o by other sets of excitations T+T+$o 
(cf Fig. 1). The second order expansion of the wave 
function tells us 



k L 



where Hi+k,i+k 



HokCi + HoiCk 

Hqq — Hi + k.i+k 



+ 



Lffl* 



H, 



i-\-k,i-\-k 



and J2\ 



(15) 



(m,n) 



runs on 



the couples T+T+$ = T^T^ $ but different from the 
couple ik. An infinite summation of diagrams would lead 



to the introduction of an EPV correction in the energy 
denominator 



EPV(i + k) = H <*Ci, 



(16) 



and to replace H 00 - H i+k , t+ k by H aa - H l+k .i+k 
EPV(i + k). 



HokCi + HoiCk + J2{m,n) (HomCn + HquC; 

Hqq - H l+k , t+ k + EPV(i + k) 



(17) 

Now taking benefit of (Eq. 14) i.e., replacing i?oi (and 
H ok ) by products Ci(H 00 - H u +EPV{i)) in Eq. 17 one 
obtains an expression closer to the CC expansion: 



Cy+^>. — C'iCk 



Hoo - H u + EPV(i) + H t 



00 



11 H - - EPV{ky\ 



+ 



E 

(m,n) 



CrnCn 



Hoo - H l+k , l+ k + EPV(i + k) 
Hoo - H mm + EPV{m) + H ao 



H nn + EPV(n) 



H 0Q - H i+k , t+ k + EPV(i + k) 



(18) 



The CC theory writes 

^T+*i = CjCfc + E] C m C ni (19) 

(m,n) 

which ignores the possible deviation of the energy de- 
nominators from additivity i.e., assumes 



Hoo - H i+ k, i+ k = Hoo - Hu + H 00 - H k k ■ Using the com- 
pact notation A, = -Hoo+Ha, A i+ k = -H o+H i+ k,i+k, 
one obtains 



C^-dCH - C t G k { - Ai+k + EPV(i + k) 1 

yr' r r f ~ Am + EPV ( m ) ' A » + EPV(n) » 
^ Cm M -A i+k + EPV(i + k) )• U0) 



(m.n) 



Notice that if there is only one route to create $ Q = <E>i+fc 
from $0 the last sum disappears. If the excitation ener- 
gies (eventually including EPV s) are additive, the full 
term vanishes in agreement with the cancellation of dis- 
connected diagrams. This will be the case in the applica- 
tions to spin lattices with nearest-neighbor interactions 
when excitations and Tj[ will concern remote bonds. 



Regarding the vectors of type2, the operators R^{^ E + ) 
which lead from $j to R^^i have to be treated pertur- 
batively as follows. Analogously to Eq. 17 one obtains 



C 



H 



,Ci 



+ , 



epv{rU) 



E 

(n,Q 



H, 



-A R+i + EPV(R?ni) 



(21) 
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C. Final equations and comments 

If P+ is a bi-electronic operator, H i R + i — H m and 
the eigenequation relative to $i will be 



(A, - EPV(i))Ci + H lQ +J2 HijCj 

, „ r ^ ( Ai-EPV^+^EPVjk) ^ f A m - EPV(m) + A n - EPVjn) s 

^ A l+/c - EPV{i + k) ^ A i+k - EPV{i + k) 

k (m,n) 

T ^*^° T+T+* =T+T+*o 

+ V + V' H m H n Ci = (} 



This will be our basic equation. It is somewhat related to 
the usual coupled cluster equations. The similarity and 
differences appear clearly if ones considers the case of an 
SCF $o- Then the vectors of the first generation are all 
the Doubles, hence the method should be compared to a 
CCD expansion. The differences are the following: 

- the effect of the Quadruples in CCD does not take 
care of possible deviations from the additivity of 
the denominators and the sum over k reduces to 

CkC\ + 2_, C m C n . (23) 

k (m,n) 
T ^*^° T+T+* =T+r+*o 

We take into account the deviation from additivity 
of energy denominators, 

- we take care of the coupling with the Singles and 
Triples through the last summation. Hence our 
method would be similar to a CCD(S,T). 

- and our (semi) perturbative treatment of the Single 
and Triples includes EPV corrections in the energy 
denominators. These corrections also appear in the 
effect of deviations from additivity of multiple ex- 
citation energies. Most of these effects were consid- 
ered in previous worksj22i2i which were expressed 
in an intermediate Hamiltonian formalism, requir- 
ing an iterative dressing of small matrices. The set 
of coupled non linear equations (Eq. 22) is much 
clearer. 



III. APPLICATIONS TO PERIODIC 1-D AND 
2-D SPIN LATTICES 

As will be shown in this section the method is of an 
extreme simplicity when applied 



- to periodic lattices of spins, (and electrons ; 22 ^ ) 

- ruled by model Hamiltonians which only introduce 
short range interactions (usually between nearest- 
neighbor atoms) 

- provided that one starts from a strongly localized 
wave function, product of atomic or bond orbitals. 

The usual implementation of the Coupled Cluster 
method are governed by a hierarchy of excitations. They 
introduce for instance all Single and Double excitations. 
Here we adopt a different strategy, since we follow the 
genealogy of the generation of from <I>o under the suc- 
cessive applications of H. Due to the localized character 
of <&o and the short range nature of the operator in H, the 
l st -generation vectors will involve short range excitations 
only, and in periodic systems one has a limited number 
of types of l st -generation excitations. Under these condi- 
tions the number of coefficients to determine is extremely 
reduced, and their values are obtained by solving a set 
of coupled polynomial equations, which may be done on 
a pocket calculator. We shall take a few examples il- 
lustrating the method and its performances, comparing 
the so-calculated values to the exact (or highly accurate) 
ones. Section A will illustrate the derivation of the SCP 
equations for simple lattices, and will compare the re- 
sulting cohesive energies to benchmarks values. Section 
B will employ the SCP method to study two phase tran- 
sition phenomena taking place in 2-D lattices. 

A. Cohesive energy of simple spin lattices 

The ruling Hamiltonian will introduce only nearest- 
neighbor interactions 

H=J2 W(SiSi - 1/4), (24) 

(id) 
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between adjacent atoms. This Hamiltonian puts to zero 
the energy of the upper multiplet 



1. 1-D chain from Neel 

If all bonds are equal, taking <I>o as the Neel function 
(spin alternation between adjacent atoms), the zero order 

energy per atom is E^ h = — J. There is only one type 



T I T I T I 



of vectors in Si, the one obtained from a spin exchange 
T- + on a bond i, of coefficient C, coupled by J with $ 
(H i0 — J) hence 



E C oh = —J(l - C), 
Ha — Hqq = Aj = 2 J, 
EPV(i) = 3CJ, 



(25) 



A simple second-order estimate of the energy would lead 



ated according to (Eq. 18) as 



_ 2 ( 2J-3CJ \ 



(27) 



= 



and the final equation is 

-3C7 2 + 2C + 1 + 2C7 2 ( IZ^c ) = °" ( 28 ) 

Its solution is C = —0.3751 hence E co h = — 1.3751J, 
which deviates by 0.8% from the exact value — 1.3862J = 
2Jln2. 



2. 2-D square lattice from Neel 



Here E. 



(o) 

coh 



-2 J, E coh = -2 J(l + C), where C is the 



coefficient relative to a spin exchange on a bond. Here 
A; = H it - H 00 = 6 J, EPV(i) = 7CJ. The coefficients 
of the doubly spin exchanged vectors are CkCi, except 
for 

- spin exchanges on the same plaquette. There are 



ti 



to C 



-1/2, E, 



(2) 
h 



1.5 J, much larger than the ex- 
Including the 
EPV' s in the denominator would give 



act value E™£ ct = -2Jln2 = -1.386J 



(2 J — 3CJ)C + J = 
3C 2 - 2C - 1 = 
C = -1/3 
E'}% = -1.333 J. 



(26) 



which is already a better estimate. There is no coupling 
between the vectors $i and <j>j € Si . The second genera- 
tion vectors are obtained by two spin exchanges on bonds 

1 and k. If the bonds i and k are far the doubly excited 
vectors T^T 4 + $o have a coefficient Ci+k = CiCk since 

- there is no other couple of l s *-generation excitation 
T+T+$ = T+T+Qo, 

- the excitation energies are additive. 

The only exceptions concern the double exchanges Tj^T^ 
and T^T^ , on second-neighbor bonds. For them Ai + k = 

2 J, EVP(i + k) = 5CJ. Their coefficient can be evalu- 



since T+T+ = T+T+ 



and 



+fe 



8J 



EPV(i + k) = Y1CJ 



hence C7 l+fe = 2C 2 (f^f) 



two such double excitations for a bond i. 

spin exchanges on other second neighbor bonds (the 
number of which is 14 for a given bond i). 



A i+k = 10J 



-h 



EPV(i + k) = 13CJ 
n — n 2 ( 12-I4c\ 

U '+ k - ° ^10-13C ) 



The final equation is 



X I i- 



T , T 



(6-7C7)C+l+2C 2 



16 - 16C 
8- 12C 



+14C 2 



2-C 
10- 13C 



= 0, 
(29) 
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the solution of which is C = -0.16327, E coh = -2.3265 J 
to be compared with the best Monte Carlo estimate 
— 2.3386J. The error is 0.51%. Our estimate may be 
compared with classical CC evaluations^^ which give 
— 2.2970J with two-body operators 
— 2.3274J with four-body operators 
— 2.3340J with six-body operators 
— 2.3364J with eigth-body operators. 
The use of many-body operators makes the algorithm 
much more complex. Our results is a modified two-body 
method, which provides the same accuracy as handling 
four-body operators. The same accuracy is obtained for 
anisotropic Hamiltonian. 



B. Identification of phase transitions 

This section illustrates the ability of the SCPE method 
to locate the critical values of physical interactions at 
which phase transition occurs. Two examples will be 
considered, both concerning 2-D lattices. In all cases a 
relevant reference function <&o is selected for each phase, 
from physical arguments. The SCP equations are estab- 
lished, leading to an evaluation of the cohesive energy for 
the corresponding phase. The step-by-step derivation of 
these equations in not given, and the equations are re- 
ported in the Appendix. The intersection of the cohesive 
energy curves as functions of the structural parameter 
enables to give a reasonable estimate of the critical value 
of the parameters. 



1. The anisotropic Heisenberg Hamiltonian for the 2-D 
square lattice 

The anisotropic Hamiltonian may be written as 

H = £ J {^S\Si + S l x Sl + SISl) . (30) 

For a 2-D square lattice and for A < — 1 the ground 
state is ferromagnetic, all sites bearing parallel spins. For 
A = — 1 the ferromagnetic state is degenerate with the 
pure XY solution, in which the bonded atoms bear al- 
ternatively (a + /3) and (a — /3) spins. This distribution 
defines a reference function $o which will be relevant for 
the — 1 < A < 1 domain. For A > 1 the wave function can 
be generated from the Neel function & in which adjacent 
atoms bear different spins. When A tends to infinity the 
system becomes Ising (^o — > $' Q ). The phase transition 
between the XY supported and the Neel supported phase 
takes place for the isotropic lattice (A = 1). Besides ana- 
lytic series expansion o 2 fii 2 7i 2 fii 2! ? accurate calculations on 
this problem have been performed according to various 
techniques, among which one may quote 

- Quantum Monte Carlo 

calculationa 3 "i 31 i 32 r S3 i 34 i 3R i 3fi i 37 i 38 



TABLE I: Cohesive energy of the anisotropic 2-D lattice. 



A 


SCP(XY) 


SCP(Neel) 


QMC 38 





-0.54672 




-0.54882 


0.1 


-0.55564 




-0.55681 


0.25 


-0.57024 




-0.57142 


0.33333 


-0.57892 




-0.57926 


0.5 


-0.59777 




-0.59832 


0.6 


-0.60965 




-0.60958 


0.75 


-0.62885 




-0.63017 


0.875 


-0.64557 




-0.64848 


0.9375 


-0.64196 




-0.65846 


0.96875 


-0.65264 




-0.66396 


1 


-0.66327 


-0.66327 


-0.66944 


1.02040 




-0.67061 


-0.67612 


1.11111 




-0.70330 


-0.70722 


1.2 




-0.73727 


-0.73920 


1.25 




-0.75696 


-0.75862 


1.5 




-0.86046 


-0.86100 


1.75 




-0.96989 


-0.96965 


2 




-1.08314 


-1.08220 



- Real Space Renormalization Group with Effective 
Interactions^ 

- Coupled Cluster expansions'^ 

- Dressed Cluster Method (DCM)^a 

The two first techniques do not introduce any reference 
function <I>o- The two last ones imply a choice of such 
reference functions either as a zero-order function for the 
Coupled Cluster Method, or as a bath in which a finite 
cluster is embedded in the DCM. Our approach belongs 
to the same family. The equations derived from $o = 
<&q Y or from $q = $N^ e l are given in the Appendix A. 
One may see that they coincide for A = 1 , which appears 
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FIG. 2: 2-D square anisotropic spin lattice, evolution of the 
cohesive energy as a function of the anisotropy parameter A. 
(o) QMC* SCP results: (...) from XY $ , (-) from Neel $(,. 
The thick symbols identify the relevant domains from $0 and 
$0 respectively. 



as the critical value of the parameter. Fig. 2 gives the 
evolution of the cohesive energy E and E' Q calculated 
from both references and a comparison with the most 
accurate results. One sees 

- that the derivatives of the energies in A = 1 are not 
identical 

(8E\- = (dEo\ , «\ = (2EY 

{dXj^ \ dx ) x=1 r { d\ J x=1 \dx)^ 

(31) 

which confirms the l s *-order character of the phase 
transition, 

- the agreement of the here-calculated cohesive ener- 
gies with the best QMC calculations in good, the 
error being lower than 1% in the — 1 < A < oo 
domain. 



2. The 1/5-depleted square lattice 

The 1/5-depleted square lattice has received interest 
in the context of the study of a real material, namely 
the CaViO^ lattice. Precise studies of that material in- 
dicated that 2™ d -neighbor interactions have to be con- 
sidered. Nevertheless the physics of the simple lattice, 
built of plaquettes and octogons, or of plaquettes con- 
nected by bonds (cf Fig. 3), appeared interesting by it- 
self. There are two types of spin interactions, namely in 
the plaquette (Jp) and in the bonds between the plaque- 
ttes (Jd)- Early studies did not give evidence of phase 
transitions, but perturbative expansions^ and Quantum 
Monte Carlo calculations^ indicate the existence of three 
phases, namely a gapped plaquette phase for J p /Jd > 
1.05 — 1.10, a gapped dimer phase for J p / Jd < 0.65 and a 
gapless Neel-type phase in the intermediate regime. We 
have reexamined this problem with the SCP equations 
starting from three different references $o> namely: 

1. a product of bond Singlets on the inter-plaquette 
bonds which is expected to be relevant in the 
Jp/Jd < 1 regime 

2. the Neel function, which should be valid in the J p ^ 
J d region 

3. a product of bond Singlets on intra-plaquette 
bonds, relevant in the J p > Jd regime. 

The corresponding equations are given in Appendix B. 
The results appear in Fig. 4 . One should remark that 
our SCP equations never diverge, even when applied in 
paradoxical regimes, for instance when one enters in the 
J P <C Jd regime from inter-plaquette bonds. The ex- 
pansion from Neel tends to an erroneous (— 2J instead 
of — 1.5J) value for Jd — (isolated plaquettes). Using 
the perturbative evaluation of the coefficient of the dou- 
bly permuted vectors (Eq. 17 instead of Eq. 18), one 
reaches the correct asymptotic value for Jd = 0. One 
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FIG. 3: The 1/5-depleted 2-D square lattice. 




Jp + Jd 



FIG. 4: 1/5-depleted 2-D square lattice. Evolution of the 
cohesive energy as a function of the J P /(J P + Jd) ratio ac- 
cording to the SCP method. (- -) from inter-plaquette bond 
Singlets "3>o, (...) from Neel <£>o, and ( — ) from plaquette Sin- 
glets "3>o- The thick symbols identify the relevant domains 
from the various reference functions. The arrows indicate the 
phase transition critical values. 



must concentrate on the estimates of the cohesive energy 
in the relevant domains. As seen from Fig. 4 one notices 
the existence of three phases. 

- a dimer phase for Jp/Jd < 0.6563 

- a Neel-order phase for 0.65 < J p /Jd < 1.0105 

- a plaquette phase for J p /Jd > 1.0105. 

Using Eq. 17 instead of Eq. 18, the critical ratios remain 
almost identical. The method confirms the existence of 
these three phases and the calculated critical values of 
the J p /Jd ratios are in good agreement with previously 
mentioned estimates. 
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IV. DISCUSSION AND CONCLUSION 

The here presented method is based on perturbative 
expansion of the wave-function but it is not strictly per- 
turbative. It is perturbative in the sense that it follows 
a genealogic generation of the wave function and that 
it uses (or refers to) pertubative estimates of the coeffi- 
cients of the second generation vectors. 
But it is not perturbative in the sense that the coefficients 
of the first generation vectors are obtained by solving a 
set of coupled polynomial equations. It may compared 
to the Coupled Cluster method (and in particular to a 
CCD(S,T) version) from which it differs by the consider- 
ation of possible deviations from additivity of multiple- 
excitation energies and by the introduction of EPV cor- 
rections in energy denominators. 

Although perfectly applicable to ab initio calculations, 
the method has been essentially designed for the study of 
periodic lattices and model Hamiltonians. In such cases, 
starting from strongly localized zero-order pictures, the 
method is analytical, it leads to simple polynomial ex- 
pressions, introducing a very limited number of variables 
(one per type of bond in most cases), which can be solved 
on a pocket calculator. The accuracy of such a simple 
method is surprising, as illustrated above. It provides an 
elegant exploratory tool for the study of spin periodic lat- 
tices. The present work has shown that the method may 



be employed to identify phase transitions. The method 
is of course applicable to Hubbard Hamiltonians, whose 
exact solutions are not known (except for the 1-D regu- 
lar chain) and to 2-D or 3-D lattices. Since the method 
only explores the physics around a given bond (up to six 
bonds in each direction) it may lead to some stimulating 
qualitative discussions to assess what is local, regional 
and properly collective in the cohesive energy of delocal- 
ized systems iSSiS^ Further work will study the spin-Peierls 
distortions. The present method essentially furnishes an 
estimate of the ground state energy. It may also provide 
short range correlation functions. The localized char- 
acter of <&o and the fact that the excitation processes 
playing a role in the low-order processes only concern a 
few bonds prevent to give estimates of the correlation be- 
yond this limit. The related Coupled Cluster formalism 
has led to specific methods for the evaluation of spectra, 
namely the Equation of Motion Coupled Cluster (EOM- 
CC) algorithm^ Future works will analyse the possibil- 
ity of such an extension in the frame of our development. 
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APPENDIX A 



SCP equations for the anisotropic 2-D square lattice (Eq. 30) 
1. equation from Neel 



(6AJ - 7JC) C + J + 2J 



/ 12AJ - UJC \ 
V 8AJ- 12 JC J ~ 



C 2 + UJ 



12AJ- UJC 
10AJ- 13 JC 



C 2 = 



E co h — J(C — — ) 



2. equation from XY 



6J-^(1 + A)C 



C+-(1 + A) + 2J 



/12J- 7J(1 + A)C\ ! 
V 8J-6J(1 + XjC J ~ 



C 2 



+14J 



12 J - 7.7 (1 + A)C 
(1 + A)C 



10 J 



13J 



c 2 - 



J 2 {l-\) 2 C 



J 2 {l~\) 2 C 



2[j-f(l + A)C] 8[J-|(1 + A)C] 



= 



(Al) 
(A2) 

(A3) 



E coh = - [(1 + A)C-1] 



(A4) 



APPENDIX B 



SCP equations for the 1/5-depleted 2-D square lattice with J p = (1 + S) and J4 = (1 — 5) 
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1. equation from inter-plaquette bond Singlets 



4(1 - 6) - !(! + <*)- ^(1 + S)C 



C+%1 + 6) 



15(1 + 5) 2 C 



4(1 - S) - 2^3(1 + 5)C 



(Bl) 



3 S A/3/. 

-2 + 2 + - (1 + <5)C 



(B2) 



2. equations from Neel function. The coefficient C (resp. C) is relative to the spin exchange in the plaquette 
(resp. between the plaquettes) 



[4 - 3(1 + S)C - 2(1 - S)C] C+ (1 + 6) + {1 + 6) 
8 - 6(1 + 6)C - 4(1 - 6)C 



t 8 - 6(1 + S)C - 4(1 - 6)C 
'4(1 - 6) - 4(1 - 6)C - 4(1 + 6)C 



- 1 



(B3) 



4(1-*) 
+2(1 + 5) 



4(1 + (5) + 2(1 - 6) - 4(1 + 6)C - 3(1 - tf)C 
4 + 4(1 + (5) - 7(1 + (5)C - 3(1 - S)C 
4(1 + 6) + 2(1 - *) - 4(1 + S)C - 3(1 - S)C 



C 2 



- 1 



CC" = 



[4(1 + S) - 4(1 + 6)C - (1 - <y)C] C + (1 - (5) 
+4(1 + J) r 4 + 4(l + J)-7(l + J)C-3(l-5)r" 



(B4) 



4(1 + ,5) 



4(1 + (5) + 2(1 - (5) - 4(1 + tf)C - 3(1 - 6)C 
%{\ + S) - %(1 + S)C- 2(1 -S)C _ 1 , 2 = 
6(1 + <5) - 6(1 + 6)C - 2(1 - 6)C 



CC' 



-\- 5 - + (l + 6)C+\(l-6)C 



(B5) 



3. equations from intra-plaquette bond Singlets. C and C concern respectively intra and inter-plaquette diexcita- 
tions 



2(1 + 6) - 2V3(1 - 6)C + (1 + 6)C] C + V3(l +6) + 



2V3(1 - 6) \±f(l - 8)C + V3(l + 8)C 



-4(1 + 5) + 2V3 [(1 - <f)C" + (1 + S)C] 



(B6) 



y (1 + -5) - ^(1 - S)C 2V3(1 + «5)C 



C + %1-6) 



(B7) 



+ - 



3(1 - 5) 2 C + 2V3(1 + <5) f^(l - 5)C + V3(l + 5)C 



-4(1 + (5) + 2V3 [(1 - 6)C> + (1 + <5)C] 



S coft = ^ -l-5+^-[(l + 5)C + (l-5)C] 



(B8) 
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